Resonant d-wave scattering in harmonic waveguides 
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We observe and analyze d-wave resonant scattering of bosons in tightly confining harmonic waveg- 
uides. It is shown that the d-wave resonance emerges in the quasi-lD regime as an imprint of a 3D 
d-wave shape resonance. A scaling relation for the position of the d-wave resonance is provided. By 
changing the trap frequency, ultracold scattering can be continuously tuned from s-wave to d-wave 
resonant behavior. The effect can be utilized for the realization of uhracold atomic gases interact- 
ing via higher partial waves and opens a novel possibility for studying strongly correlated atomic 
systems beyond s-wave physics. 
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I. INTRODUCTION 

The control of the interaction via resonances in atomic 
and molecular many-body systems is of major importance 
in the field of ultracold quantum gases. Magnetic and 
optical Feshbach resonances are well estabhshed tools of 
experimentally tuning the interactions of atomic ensem- 
bles recent years, however, an alternative way 
of controlling the inter-particle correlations has attracted 
increasing interest: the confinement-induced resonance 
(CIR) [J-[y]. It emerges due to the interplay of atom- 
atom scattering with the tightly confining geometry of 
the trap, and yields drastic changes of the scattering pro- 
cess and its properties 0, Q . CIRs have been extensively 
studied e.g. in the context of three-body [H, [9| and four- 
body scattering in a confining potential, of p-wave 
scatte ring of fermions [ll| , of scattering in mixed dimen- 
sions [l2[ or multichannel scattering in cylindrical con- 
finement [l^. For distinguishable atoms, two major ef- 
fects are the so-called dual CIR yielding a complete sup- 
pression of quantum scatterin g, a nd resonant molecule 
formation in tight waveguides |16l [l7| . Remarkable ex- 
perimental progress has lead to the observation of CIRs 
for both bosons [6|, El, [3 fermions [5] , relying on the 
control of s-wave and p-wave interactions in bosonic and 
fermionic low-dimensional systems, respectively. It is an 
intriguing perspective to extend the resonant scattering 
physics to higher partial wave interactions in quantum 
gases, since the latter is expected to provide novel many- 
body phenomena. A promising example is the possibil- 
ity that higher partial wave interactions are responsible 
fo^h^h-temperature superconductivity and superfluidity 

In the present work we suggest a novel mechanism for 
resonant higher partial-wave interatomic interactions in 
bosonic quantum gases, based on the effect of resonant 
d-wave interactions for ultracold collisions in tightly con- 
fining waveguides. We provide a resonance condition im- 
plying the dependence of the position of the d- wave res- 
onance (DWR) on the oscillator length a± = ^/h/ {iiuoi_) 



{uj±_ is the confinement frequency) of the harmonic waveg- 
uide. It is shown that this condition does not depend on 
the short-range part of the interatomic interaction but on 
the long-range part via the dispersion coefficient Cq . Ad- 
ditionally, the strength of the d-wave interactions can be 
controlled by changing the width of the waveguide, as 
it is the case for s-wave scattering. Due to the fact that 
DWRs occur in the background of s-wave CIRs, an ad- 
justment of the trap width can change entirely the nature 
of the underlying interactions of the colliding atoms from 
s- to d-wave and vice versa. We find that the DWR ob- 
served here is based on a shape resonance formed in the 
presence of both the centrifugal barrier and the confining 
trap. 



II. HAMILTONIAN, METHODOLOGY 
AND SETUP 



We consider ultracold collisions of identical bosons in a 
harmonic waveguide with a transverse potential ^iiuj\p^ 
{p = r sin 0). The latter permits a separation of the center 
of mass and relative motion yielding the following Hamil- 
tonian for the relative motion: 



2(1 r2 



2 , vir) (1) 



V{r) is the interatomic potential, where r is the relative 
radial coordinate and /i = m/2 is the reduced mass of 
the two bosons. The boundary conditions for quasi- ID 
scattering in the waveguide read, for z = rcosO ±oo. 
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$0,0 (r sin 0) (2) 
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where fo is the elastic scattering amplitude of the 
ground transversal channel, <I>o,o(^ sin 6>) is the ground- 
state wave-function of the 2D harm onic oscillato r with 
transversal energy e± = JiijJ±^ k = ^/2jj.{e — Jiw±) is the 
relative momentum of the colliding pair, and e = s\\ -\- s± 
is the energy of the relative two-body motion written as 
a sum of the longitudinal e\\ and transversal e± energies. 
We focus on the single mode regime for which the en- 
ergy lies below the first excited transversal energy level, 
Jiijj± < £ < 3Jiw±. Obviously, the scattering state ([2|) is 
parity (r — r) symmetric, corresponding to the case of 
two colhding identical bosons. 
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FIG. 1: (color online) (a) Partial cross sections for s- and d- 
wave scattering in free space, (b) transmission coefficient T in 
the waveguide and (c) interaction strength gm as a function 
of the parameter C12 (Ce = 28.82 for Cs). 



The Hamiltonian allows for a separation of the az- 
imuthal motion, corresponding to the conservation of 
the projection of the angular momentum onto the sym- 
metry axis of the waveguide. Consequently, we obtain 
a 2D scattering problem [Eqs. ([T]) and(|2])] in spherical 
coordinates, which is solved by em ploy ing the discrete- 
variable method suggested in Refs. [Isl, |23|, with the ra- 
dial part of the Schrodinger equation discretized using 
a B-spline basis [3, [HI. The interatomic interaction 
V{r) is modeled via the Lennard- Jones 6-12 potential 
V{r) = Ci2/r^^ — Cq/v^ . We hereafter employ the units 
mcs/2 = = cjo = 1, where mcs is the mass of the Cs 
atom and ujq = 27r x 10 MHz in SI units. The following in- 
vestigations are performed for the dispersion coefficients 
Cq defining the long-range part of the interatomic po- 
tential for three different atomic species, Cs, Rb and Sr 
[2i,[29|. The longitudinal energy is set to e\\ = 2 x 10"^ 
and the transversal energy is varied within the interval 
2 X 10~^ ^ £± < 2 X 10~^, corresponding to a range 
47r KHz < uj± < 4007r KHz (SI) for the waveguide con- 
finement frequency. We thereby focus on the low energy 
regime, characterized by ka± <C 1. 

Important quantities of our analysis are the scattering 
amplitude /o and the transmission coefficient 

r=|i + /o|' (3) 

which we analyze in the limit of very small longitudinal 
collision energies for which the s-wave CIR was initially 
defined 
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FIG. 2: (color online) (a) Partial cross section a for d-wave 
scattering in free space and (b) transmission coefficient for 
several values of the confinement frequency (jJ±, as a function 
of the parameter C12 around the d-wave resonance. 



III. RESULTS, ANALYSIS AND DISCUSSION 

Let us first briefly address free space Cs-Cs collisions 
(Ce = 28.82 [28|,[29|) in the energy range for which s-wave 
scattering dominates the elastic scattering process. By 
varying C12 (which changes the short range part of the 
interatomic interaction) we observe s- and d-wave res- 
onances in the corresponding partial cross sections [see 
Fig. [TJa)]. As expected, the width of the d-wave res- 
onance is by orders of magnitude smaller than that of 
the s-wave resonance. In the presence of the waveg- 
uide [see Fig. H^b)] the free near resonant s-wave scat- 
tering turns into a s-wave confinement-induced resonance 
0, 05 [m, [23, [HI, which appears as a broad dip in 
the transmission coefficient with a zero- valued minimum 
at C12 = 1.65, coinciding with a divergence of the inter- 
action strength guj = lim^u ^0(^1^6 /o/Im/o). In order 
to compare the numerically calculated position of the s- 
wave CIR with the analytical result Q we map the C12 
parameter onto the s-wave scattering length in free 
space and obtain the ratio as{Ci2 = 1.65) /a± = 0.65 
which is in good agreement with the analytical predic- 
tion as/a± = 0.68. 

The d-wave resonance persists also in the waveguide, 
where it appears as a strongly varying transmission coef- 
ficient between the limiting values and 1 on top of the 
s-wave scattering background [Fig. [Hb)]. As in the case 
of the s-wave CIR, we observe a divergence of giD at the 
resonant value C12 = 0.9988731 [Fig. [U^c)]. 

The transversal confinement leads to a shift of the d- 
wave resonance when passing from free space to confined 
scattering, as is clearly seen in the high resolution graph 
of Fig. [2j This shift, as well as the corresponding width, 
increase with the trap frequency uj± [see Fig. [2jb)]. We 
also observe a strong suppression of the s-wave back- 
ground with increasing ujj_. The above behavior of the 
transmission coefficient T{uj ±^€12) under the action of 
the confining potential can be interpreted in terms of a 
strong coupling of s- and d- waves in harmonic traps, since 
the trap potential can be represented as a sum of the Leg- 



3 



20 
10 

-10 
-20 



■ 


ad \ 

^ 


■ 


\ 



0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 

FIG. 3: (color online) s- and d-wave scattering lengths as and 
ad in free space as a function of the C12 parameter, where 
ad is defined through af+^ = - tan [1^ for ^ = 2. 

ai^^ {uj±_ = 0) = a^^Q is the value of the s-wave scattering length 
at the point of divergence of the d-wave scattering length ad- 



endre polynomials: 



i/ic^ip' = i/ia;ir2[2Po(cos^) -P2(cos^)] 



(4) 



Consequently, the confining potential induces a coupling 
between partial waves ^ and ^ + 2 in the course of the 
atomic collisions in the trap. This coupling can be tuned 
by changing the trap frequency. 

We now analyze more precisely the condition for the 
occurrence of the DWR. The experimental determination 
of the d-wave scattering length and in particular of the 
position of the corresponding resonance is a difficult task. 
However, in the following we will show the existence of a 
very useful scaling relation for the position of the DWR 
that provides us with an accurate estimate of its position. 

Let us firstly introduce the quantity ai^\ which is the 
free space s-wave scattering length at the position of 
the free space d-wave shape resonance. This is well de- 
fined since in general varies smoothly in the vicin- 
ity of the d-wave resonance (see Fig. [3j). Similarly in 
[23 the crucial role of for the analysis of the ultra- 
cold scattering near the shape resonances in free space 
was shown. We denote its value a^f^iuj^ = 0) in free 
The next step is to determine the shift 



space by q . 



Aa, 



Ad) 



of the position of the DWR in the 



waveguide a^f^ = as'^\uj± ^ 0) with respect to the reso- 
nance position in free space, in the limit e\\ 0. Note 
that as is altered here by varying the C12 parameter, 
while experimentally it may be achieved by employing 
e.g. magnetic Feshbach resonances [3, 6]. In Fig. SJ^a) we 
show the dependence of the quotient p = Aai^^ /Cq on 
for the three different atomic s pec ies Cs, Rb and 
Sr with corresponding Cq coefficients [281 , |29[ . By fitting 
the numerical data we find this dependence to be rather 
accurately described by the linear relation 
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for a comparatively broad range of values of uj± for all 
considered Cq coefficients. This shows that the quotient 
p changes with the confinement strength independently 
of the type of colhding atomic species. In contrast to the 
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FIG. 4: (color online) Quotient p = Aai^^ /Cq as a function of 
l/a\ (a) for different values of Cq corresponding to the atomic 
species Cs, Rb and Sr and (b) for Rb atoms with different 
interatomic potentials (see text), for confinement frequencies 
47rKHz < uj± < 4007rKHz. Dashed lines correspond to the 
resonance condition given in Eq. ((S]). 



s-wave CIR, whose position depends exclusively on the ra- 
tio as/a±^ the shift Aai^"^ of the DWR depends addition- 
ally on the dispersion coefficient Cq. This is because the 
Van-der-Waals potential tail determines the width of the 
centrifugal barrier (absent in s-wave scattering), which 
in turn strongly influences the DWR state. Note that p 
possesses the correct asymptotic value p{a± 00) = 
in the limiting case of scattering in free space. Addition- 
ally, a fitting of experimental data with the Eq. (j5j) can 
provide an alternative way for measuring the dispersion 
coefficient Cq. 

In order to assure that Eq. ([5]) is not specific to the 
model used for the interatomic interaction, we have fur- 
ther studied a hybrid potential which possesses a short- 
range part different from the Lennard-Jones 6-12 poten- 
tial, reading: 



Vh{r)=Ae- 



■qr 



(6) 



where A and q are constants and fcix) is a cut-off function 



fc{r) = 9{r - Tc) + 9{rc - r) exp 



(7) 



with 0{x) being the Heaviside step function and Vc the 
cut-off radius. We choose A = 5000 and Tc = 1.2 and vary 
the parameter q. Fig.Hl^b) shows p for Rb-Rb collisions as 
a function of 1/a^ for the Lennard-Jones and the hybrid 
potential. The values for p{l/a\) are for both potentials 
in very good agreement with the scaling relation Eq. (|5]). 
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FIG. 5: (color online) Probability density \^{r,0)\^ for Cs-Cs 
collisions (a) for the off-resonant case, (b) at the s-wave CIR 
and (c) at the DWR, for cj^ = 2 x 10"^. 
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FIG. 7: (color online) Probability density \^(r,0)\^ for Cs- 
Cs collisions (a) a; = 5.7 x 10"^ ('left side' of the DWR), 
(b) uj = 5.96 X 10"^ (position of the DWR) and (c) uj = 
6.36 X 10-^('right side' of the DWR). 



Let us now inspect the probability density |^(r, ^)p 
for different regimes of Cs-Cs collisions, plotted in Fig. [H 
The non-resonant case [Fig. [5ja)] leads to a probability 
density which is substantial only far from the origin r = 
0. Fig. ^h) and (c) refer to the s-wave CIR and DWR 
in the waveguide, respectively. For the s-wave CIR we 
observe a localization of the probability density near the 
origin, while the angular distribution possess no nodes. 
In contrast, the DWR exhibits, as expected, two nodes in 
the coordinate. 

Fig. [6] together with Fig. [7] illustrate how the d-wave 
interactions in harmonic waveguides can be controlled by 
altering the trap width. Varying uo±_ (a^) changes the 
transmission coefficient T near the DWR and accordingly, 
and so the strength of the d-wave interaction between the 
bosons can be tuned. In order to further illustrate this we 
show in Fig.[7fa)-(c) the probability density distributions 
for the resonant value of uj^ as well as far from the reso- 
nant region. A weak anisotropy is observed far from reso- 
nance whereas strong anisotropy characterizes the density 
profile close to resonance [Fig. [Tl^b)]. We emphasize the 
principal difference in the resonant behavior of the trans- 
mission coefficient near the s-wave CIR and DWR [see 
Figs.[Hb) and [6]. Near the s-wave CIR the contribution 
of d-wave scattering is negligible with respect to the total 
ID scattering amplitude. However, near the DWR the 
d- and s-wave scattering states are strongly coupled and 
their interference lead to a sharp variation in the trans- 
mission coefficient between unity and zero, characterized 
by the Fano asymmetric lineshape [31]. Changing the 
confinement length one can therefore completely alter 
the nature of the bosonic interactions from s- to d-wave 
character and vice versa, since the s- and d-wave scatter- 
ing coexist in the transmission coefficient. From Fig. [71 



one can conclude that the region of significant d-wave 
scattering is rather broad with respect to the variation 
of uj±_ and exceeds the width of the d-wave resonance in 
Fig.El 

To clarify the physical reason for the appearance of 
the DWR let us analyze the nature of the free space and 
confinement-resonant states. Fig. [2fb) provides us with 
the values of the C12 parameter for each uo^ at resonance, 
which we then use in the 3D scattering problem to extract 
the corresponding free space resonant energies E^^\ It 
is found that these coincide with the the ground state 
energy in the harmonic confinement with increasing C12'. 
Er^^^ = huj± (see Fig. [8]). This allows us to the conclude 
upon the stability of the position of the near-threshold d- 
wave resonance Er^^"^ while adding the trap potential (see 
Fig. [9] for the following discussion). Behind the centrifu- 
gal barrier (r < tq) the trapping potential only slightly 

shifts the interatomic spectrum Er^^^ Er^^^ ~ Er^^^ 
for ro <C a^. However, in front of the centrifugal barrier 
(r > ro), where V{r) 0, the trap potential becomes 
dominant and leads to the quantization of the interatomic 
continuum. 

The near-threshold d-wave resonance Er^^^ in free 
space {uj± = 0) becomes a weakly-bound state Er^^"^ ^ 
with binding energy |£^r^^^ — huj±\ ~ huj± referring as a 
threshold to the first vibrational state in the waveguide. 
Thus, it becomes clear that for transforming this weakly- 
bound state into the DWR it has to be shifted in order to 
lead to an energetically degeneracy with the threshold of 
the transversal ground state, yielding then Er^^^ — huj±. 
The latter is here achieved [see Fig. [2jb)] by altering the 
C12 coefficient. Experimentally, it can be realized by the 
magnetic Feshbach resonance technique. Further, one can 
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FIG. 6: (color online) Transmission coefficient for C12 = 
0.9989621 as a function of the confinement frequency u± near 
the DWR. 
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FIG. 8: (color online) Ratio of the energy Er^^^ of the res- 
onant state in free space and the confining frequency (jJ± as 
a function of the C12 parameter. u± is varied in the range 
47r KHz <0J±< 4007V KHz. 
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FIG. 9: (color online) A schematic illustration of the scat- 
tering potential in free space (blue) and in the presence of 
confinement (black dashed). 



properties in detail. The DWR observed here is a shape- 
resonance as opposed to the Feshbach resonance of the 
s-wave CIR. It is controllable by adjusting the frequency 
of the trap. The confinement plays here a two- fold role: 
It can control the strength of the interactions among the 
bosons as well as alter the nature of the interactions from 
s- to d-wave character. The DWR effect can be analyzed 
and exploited in experiments analogous to [ssf with the 
help of the scaling relation (|5]). The observation and ap- 
plication of it in anisotropic waveguides represents an in- 
triguing perspective, since in the limit of the quasi-2D 
regime the anisotropic profile of the d-wave interactions 
can be exploited to introduce novel properties into quan- 
tum gases, or more specifically speaking, strongly corre- 
lated bosonic systems. Unconventional superfluidity and 
superconductivity in the presence of higher partial wave 
interactions are examples supposed to bear important pe- 
culiarities. 



control the appearance of the DWR exclusively by tuning 
the trap frequency uo±_ to arrive at the resonance condition 
E^^^ Er^^^ = huo± without altering the interatomic 
interaction V{r). The above analysis also explains the 
increase of the width of the DWR with increasing C12 [see 
Fig.[2{b)] due to a corresponding decrease of the relevant 
width of the centrifugal barrier. Finally, it demonstrates 
that, in contrast to the s-wave Feshbach CIR [32| which 
arises from a bound state of the first excited transversal 
channel, the observed DWR is a shape resonance. 

IV. BRIEF SUMMARY 

We have demonstrated how strong d-wave interactions 
can develop in a harmonic waveguide and analyzed their 
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